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' Abstract. We performed an observational program with the X-ray satellite BeppoSAX , to study objects with extreme syn- 

chrotron peak frequencies (i^pcak > 1 keV). Of the seven sources observed, five revealed or confirmed their extreme na- 
ture. Four showed peak frequencies in the range 1-5 keV, while one (lES 1426-1-428) displayed a flat power law spectrum 
(Ox = 0.92 ± 0.04) which locates its synchrotron peak at or above 100 keV. This is the third source of this type ever found, 
^ ^, after Mkn 501 and lES 2344-1-514. In the context of the whole blazar class, the broad band properties of these objects confirm 

I ' the scenario of a synchrotron peak smoothly spanning the IR - X-ray range, which explains the multi-frequency properties of 

O I the blazar class. Our data also confirm the large i^pcak variability which seems to characterize this class of sources, compared 

■ with lower fpcak objects. Given the high synchrotron peak energies, which flag the presence of high energy — ^ — " fu„„„ 

^ ' extreme BL Lacs are also good candidates for TeV emission, and therefore good probes of the IR background. 



Key words. BL Lacertae objects: individual: lES 0033-1-595, lES 0120-1-340, PKS 0548-322, GB 1114-F203, lES 1218-1-304, 
lES 1426-1-428, H 2356-309 - X-rays: general - TeV: general 

1 . Introduction Padovani & Giommi 1995). In the past the picture given by X- 

ray and radio surveys presented a bimodality in the distribution 



Among active galactic nuclei (AGNs), blazars are the objects 
most dominated by non thermal continuum emission, which ex- 
tends from radio to gamma rays, and whose properties are best 
explained by emitting plasma in relativistic motion towards the 
observer, closely aligned with the line of sight (see e.g. UiTy & 



of the sources, with radio-selected objects mostly LBLs, and 

X-ray selected mostly HBLs. New results are now revealing a 

more continuous scenario. On one hand, new surveys like the 

DXRBS (Perlman et al. 1998), the RGB (Laurent-Muehleisen 

T> A • ioo<\ T-u • .1 A- . -u /ccT^^ • ct al. 1998) and thc REX (Caccianiga et al. 1999) afc showlng 

Padovani 1995). Their spectral energy distribution (SED), in a . j ^ i „ • i 

^ ... • , , • , , , a more continuous and smoother range in peak rrequencies. 

i^i'i/ representation, IS characterized by two main broad peaks: „, „■ , nr-i^ • • ^„ j , , .1 

, f . , , . , , . The gap in the SEDs properties is now nlled thanks to the 

that at low energy IS commonly explained as synchrotron emis- ■ r ■ , 1 j • ^ 1 

, . , , , . . . , ^ , ., ; , . , covering or previously unexplored regions or the parameters 
sion by highly relativistic electrons, while that at higher en- /•t^t^ji.ji.j j 

■' o space (in a;, i'r and broad band spectral indices ttrx, ttr o and 

ergy is probably due to Inverse Compton emission of the same x u j r. o a i. r » ji 

, , , , . , , ., , Qfox)- On the other hand, the BennoSAX observations or Mkn 

electrons off seed photons, which can be either the same syn- cm /d- , 1 inm\ ^ ice oa/i/i , • » 1 oAr>A\ 

, „ , r , 501 (Plan et al. 1997) and lES 2344H-5 14 (Giommi et al. 2000) 

chrotron photons produced locally, or photons from the exter- , 1 j . u . 1.1 u 

^ ^ ±- have revealed that the synchrotron peak not only can reach the 

nal environment (see e.g. Ghisellini et al. 1996; Sikora et al. ^ , „ „ ... 

2-10 keV band (as was known tor few sources, like Mkn 421), 

1997; Ghisellini &Madau 1997; DermeretaL 1997). , , I ^, , mn. w . ■ « ■ 

but can go even further, up or above 1 00 ke V, at least in a flaring 

Among blazars, BL Lacertae objects are the sources with . . rp, 1 u u . j j u 1 

° ■' state. These exceptional observations have extended by nearly 

the highest variety of synchrotron peak frequencies, ranging , . ^ j 1 " 1 ui " f 

u Tr. 1 u iT\r f V u J / 11 J T two ofdcrs of magnitude the known available range of syn- 



from the IR-optical to the UV-soft-X bands (called Low or 
High energy peak BL Lacs, i.e. LBL or HBL, respectively; see 



chrotron peak frequencies, disclosing a new region of possible 
physical parameters that can give us valuable insights on the 
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At present, however, very little is known about such ex- 
treme objects: we do not know if they are very rare or not, 
and if these energies are reached only during flares or also in 
normal, "quiescent" states. Another interesting topic is if there 
exists an upper limit to the peak frequency, and if we have al- 
ready observed it or not. In order to answer to these questions, 
a necessary step is to find and study other objects of this kind, 
sampling more accurately the high energy region (> 1 keV) of 
the peak sequence. With this goal,we have proposed and per- 
formed an observational campaign using the BeppoSAX satel- 
lite (Boella et al. 1997), after selecting a sample of candidates. 
The X-ray band is ideal for this purpose: in this band, objects 
characterized by different synchrotron peak frequencies show 
different spectral properties, according to the dominant radia- 
tion process (synchrotron or inverse Compton). This provides 
means to roughly estimate the position of the synchrotron peak, 
with the aid of the overall SED shape to recognize the origin of 
the X-ray emission. 

In fact, the X-ray spectral index might even be considered 
the parameter by which to classify the source: BL Lacs which 
have the synchrotron peak in the IR band (i.e. LBLs), have a 
flat (ax < 1) X-ray spectral index, due to inverse Compton 
emission, as for FSRQs []. As the peak shifts towards higher 
energies, the synchrotron emission becomes more and more 
dominant in the X-ray band, and so we have, in sequence: in- 
termediate BL Lacs, where we see the high energy tail of the 
synchrotron emission in the soft X-ray band and the flat inverse 
Compton spectrum in the hard band, resulting in a concave X- 
ray spectral shape (agoft > 1 and ahard < 1); HBLs, where 
the X-ray band is completely dominated by the synchrotron 
emission above its peak, and so show a steep spectral index 
(ax > 1); and Extreme BL Lacs, where the synchrotron peak 
becomes visible directly in the X-ray band at energies > 1 keV 
(resulting in a convex spectrum with cksoft < 1 and ahard > 1, 
or in a flat X-ray slope all the way up to the maximum ob- 
servable X-ray energies, as in the case of Mkn 501). Thanks 
to its wide energy band (0. 1-100 keV), the BeppoSAX satellite 
offers the best opportunity to identify and study these extreme 
sources. 



2. The sample 

Our candidates have been selected from the HEAOl survey, 
the Einstein Slew Survey and the Rosat All Sky Survey Bright 
Sources Catalogue (RASSBSC), using also information in 
other bands to construct the SEDs and to determine the broad 
band spectral indices (a^o, cnox and arx)- These are good in- 
dicators of the location of the synchrotron peak, since objects 
of different characteristics tend to gather in different regions of 
the parameter space (Stocke et al. 1991; Padovani & Giommi 
1995; Fossati et al. 1998; Padovani et al. 1997; see also Fig. 1). 

The selection criteria were based on properties suggesting 
a high value of the synchrotron peak energy: 
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' We call flat or steep a power law spectrum with index a < 1 or 
a > 1 respectively. This convention corresponds to rising or decaying 
(with frequency) spectra in a u-vF^ plot. 



Fig. 1. The plane aox-c^rx for the SLEW sample of BL Lacs, 
the 1 Jy sample of BL Lacs and the 2 Jy sample of FSRQs. All 
sources with arx < 0.7 are X-ray selected BL Lacs. The SED 
of all sources has been fitted with a cubic function in order to 
find the synchrotron peak frequency v-pcak (Fossati et al. 1998). 
Filled squares mark the location of the observed candidates. 



- a very high Fj-Zi^iadio ratio (> 3 x lO^^*^ erg cm^^ s^^/Jy, 
measured at least once, in the [0.1-2.4] keV band and at 5 
GHz respectively); 

- flat X-ray spectrum (when available), connecting smoothly 
with the flux at lower frequencies; 

- appropriate values of a^o, ctox and arx (see e.g. Fig. 1). 

In addition, a high X-ray flux (> 10~^^ erg s^^ cm^'^) was 
also requested, to achieve a good detection in the PDS instru- 
ment. 



3. BeppoSAX observations 

In Table 1 we report the observation log and some basic in- 
formation for the seven observed sources. The LECS and PDS 
exposure times are usually shorter than the MFCS ones, due to 
light leakage in the LECS instrument, which forces to take the 
exposures only during the part of the orbit in the earth shadow 
cone; and due to the on-off source switching for the PDS in- 
strument. 

The data analysis was based on the linearized, cleaned 
event files obtained from the BeppoSAX Science Data Center 
(SDC) on-line archive (Giommi & Fiore 1997). The data from 
the single MFCS units were merged in one event file by SDC, 
based on sky coordinates. The event file was then screened with 
a time filter (from SDC) to exclude those intervals related to 
events without attitude solution (i.e. conversion from detector 
to sky coordinates), as recommended by SDC (see Handbook 
for NFI Spectral Analysis, Fiore et al., 1999). Due to remaining 
calibration uncertainties, the channels 1-10 and over 4 keV for 
the LECS, and 0-36 and 220-256 for the MFCS were excluded 
from the spectral analysis. Standard extraction radii of 4 and 8 
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Table 1. Log of the observations 



lES 1426+428 



bintime 10000 s 



Date 



Detector Exp. Time Net countrate full band 
seconds counts/s 



lES 0033+595 


z= 


9 


gal. Nh= 


4.24 xlO^^ cm 


2 


18-19 Dec. 99 


LEGS 
MECS 
PDS 


13267 
43424 
22877 




0.368 ±0.006 
0.766 ±0.004 
0.503 ±0.052 




lES 0120+340 


z= 


0.272 




= 5.15 xlO^° cm 


-2 


3^ Jan. 99 


LEGS 
MEGS 
PDS 


9661 

32039 

16263 




0.198 ±0.005 
0.260 ±0.003 
0.381 ±0.065 




2-3 Feb. 99 


LEGS 
MEGS 
PDS 


4848 

21562 

10128 




0.179 ±0.008 
0.213 ±0.003 
0.196 ±0.082 




PKS 0548-322 


z= 


0.069 


gal. Nh= 


= 2.51 xlO^° cm 


-2 


20 Feb. 99 


LEGS 

MEGS 
PDS 


9328 

12439 

9355 




0.280 ±0.006 
0.350 ±0.005 
0.142 ±0.083 




26-27 Feb. 99 


MEGS 


2025 




0.309 ±0.013 




7-8 April 99 


LEGS 
MEGS 
PDS 


5251 

18943 

10164 




0.222 ±0.008 
0.248 ±0.004 
0.039 ±0.085 




GB 1114+203 


z= 




gal. Nh= 


1.36 xlO^° cm" 


2 


13-15 Dec. 99 


LEGS 
MEGS 
PDS 


13551 
42355 
20432 




0.195 ±0.004 
0.128 ±0.002 
-0.026 ±0.051 




lES 1218+304 


z= 


0.182 


gal. Nh= 


1.78 xlO^° cm 


-2 


12-14 July 99 


LEGS 
MEGS 
PDS 


10626 

42693 
20670 




0.284 ±0.006 

0.268 ±0.003 
0.108 ±0.051 




lES 1426+428 


/.= 


0.129 


gal. Nh= 


1.36 xlO^° cm" 


-2 


8-9 Feb. 99 


LEGS 
MEGS 
PDS 


15153 

40657 
20432 




0.232 ±0.004 

0.293 ±0.003 
0.487 ±0.058 




H 2356-309 


z= 


0.165 


gal. Nh= 


= 1.33 xlO^° cm 


-2 


21-22 June 98 


LEGS 
MEGS 
PDS 


15098 
41034 
18035 




0.327 ±0.005 
0.374 ±0.003 
0.150 ±0.063 





arcmin for MECS and LECS respectively were used, except for 
PKS 0548-322 (see paragraph 3.2.3). 

The PDS instrument has no imaging capability, and its 
f.o.v. (~ 1.4° FWHM) is larger than the other Narrow Field 
Instruments (radius ~ 28' for the MECS). Therefore there is 
the possibility for PDS spectra to be contaminated by hard 
serendipitous sources in the f.o.v, not visible in the smaller 
MECS images. To explore this possibility, for each object we 
checked in the NED and WGACAT databases for the presence 
of potentially contaminating sources. 




10 15 
Time (hours) 



Fig. 2. MECS lightcurve and Hardness Ratio for lES 
1426+428, defined as (5-10 keV)/(2-5 keV) countrates. 



3.1. Time analysis 

We looked for time variability in each observation, using time 
bins from 500 to 10,000 s: no significant variability has been 
detected (x'^-test Constance probabilities > 20%), except for 
lES 1426+428. For this source the lightcurve clearly shows 
a trend (see Fig. 2), with a flux decreasing in the central part 
of the observation by about ~ 10% in 10 hours, before rising 
again. This variabiUty is indicated by the low x^-test probabil- 
ity of constant flux (< 0.02%), and confirmed by the fact that 
the same pattern is visible in the lightcurves from the single 
MECS2 and MECS3 instruments. No significant variations are 
present in the LECS lightcurves, but the lower statistics does 
not allow to draw conclusions on a different behaviour of the 
low energy emission. Fig. 2 also shows the corresponding hard- 
ness ratio (HR): the scatter is not compatible with a constant 
value. After the 1 1th hour there seems to be an anticorrelation 
between the HR and the countrate. However the significance of 
this anticorrelation is small: there is still a ^20% probability to 
obtain the same correlation coefficient from a random sample. 

3.2. Spectral Analysis 

The spectral analysis was performed with the XSPEC 10.0 
package, using the latest available response matrices and 
blank-sky background files (at 01/2000). 

Particular attention was paid to the LECS blank-sky files, 
chosing those extracted in the same coordinate frame as the 
source files, to avoid the errors in the background subtraction 
arising from the different scaling of the area between raw and 
sky-detector coordinates. 

The spectra were rebinned both with more than 20 net 
counts in each new bin and using the rebinning files provided 
by SDC, applying the Gehrels statistical weight (Gehrels 1986) 
when the resulting net counts were below 20 (typically 12-15 
in some low energy channels of LECS). Various checks have 
shown that our results are independent, within the uncertain- 
ties, of the adopted rebinning. The LECS/MECS normalization 
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Single power-law fits, Lecs+Mecs 



Source 


Date 


Nh 




-Fl keV 


F2- 


10 keV 




Lees 


X?/d.o.f. 


notes 


(lAU name) 




10'" cm-'' 




/xJy 


ergs cm 


s-i 


Mecs 






0033+595 


18/12/99 


42.4 fixed 
75 ±7 


0.93 ± 0.02 
1.05 ± 0.04 


13.5 ±0.3 
16.8 ± 0.7 


5.9 
6.0 


xlO" 
xlO" 


-11 
-11 


0.60 
0.69 


2.24/152 
1.12/151 


±PDS 


0120+340 


3/1/99 


5.15 fixed 

8-5+4.3 

Q Q-1-9 
^•^+3.0 


1.08 ± 0.04 
1.11 ±0.06 
i.io ± U.Ud 


5.0 ±0.3 
5.2 ± 0.3 

D.Z ± U.O 


1.7 
1.7 
1.7 


xlO" 
xlO" 
xlO" 


-11 
-11 
-11 


0.68 

0.70 
n '71 


0.98/122 
0.96/121 
u.yo/izz 


LECS> 0.3 keV^' 
LECS> 0.3 keV^' 
LxiL-o total ^ 




2/2/99 


5.15 fixed 

S Q-2-4 
°-y+6.1 


1.23 ± 0.05 
i.o± ^ u.uy 


4.8 ±0.3 

5-4+0.5 


1.3 
1.3 


xlO" 
xlO" 


-11 
-11 


0.63 

fi Aft 
U.Oo 


1.15/95 
n OA /OA 




0548-322 


20/2/99 


2.51 fixed 

A Q-0-8 

4-y+i.i 


0. 88 ± 0.04 

1 no _L n fiV 

1. Uo ± U.U ( 


5.1 ±0.3 

D.o ± U.4 


2.4 
2.3 


xlO" 
xlO" 


-11 
-11 


0.77 

\j. 1 J 


1.83/85 






26/2/99 


2.51 fixed 


1.32 ±0.19 


7 5-1-5 

'•3+1.9 


1.8 


xlO" 


-11 




0.58/17 


MECS only 




7/4/99 


2.51 fixed 

A 7-10 
^•'+1.6 


1.14 ±0.05 
1.24 ± 0.08 


5.0 ±0.3 
5.7 ±0.4 


1.6 
1.5 


xlO" 
xlO" 


-11 

-11 


0.66 
0.68 


1.10/79 
0.72/78 




1114+203 


13/12/99 


1.36 fixed 
3.1 ±0.4 


1.53 ±0.03 

-■-•8^+0.05 


3.7 ±0.2 

^•^+0.2 


6.7 
6.4 


xlO" 
xlO" 


-12 
-12 


0.74 
0.65 


3.14/83 
1.42/82 




1218+304 


12/7/99 


1.78 fixed 

o ri-O.S 
"'•"-'+0.4 


1.39 ± 0.03 
1.51 ±0.04 


7 1-0.3 
'•1+0.2 

8.1 ±0.3 


1.6 
1.5 


xlO" 
xlO" 


-11 
-11 


0.66 
0.67 


1.95/94 
1.15/93 




1426+428 


8/2/99 


1.36 fixed 

^■•-'+0.4 


0.91 ±0.03 
0.92 ± 0.04 


4.6 ± 0.2 
4.6 ±0.2 


2.0 
2.0 


xlO" 
xlO" 


-11 
-11 


0.68 
0.68 


LOO/90 
LOO/89 




2356-309 


21/6/99 


1.33 /ixed 
2.3 ± 0.3 


0.96 ± 0.02 
1.06 ± 0.04 


6.1 ±0.2 
7.0 ±0.3 


2.5 
2.5 


xlO" 
xlO" 


-11 
-11 


0.72 
0.70 


3.17/37 
1.58/36 


+ PDS 



Errors at 90% conf. level for 1 (fixed Nn) and 2 (free .^Vh) par. of interest. see text 



Table 2. Single power-law fits with galactic and free absorption. In bold are highlighted the dates of those observations for which 
the single power-law model provides an adequate representation of the data. For the other observations, instead, a significantly 
better fit is provided by the broken power-law models (see Table 3). 



factor was let free to vary: normal acceptable values should 
be within the 0.65-1.0 range, as indicated by SDC. For the 
PDS instead, the normalization factor was fixed at 0.9, as re- 
sults from intercalibration tests performed on known sources 
(Fiore et al. 1999). 

In order to unveil the possible presence of the synchrotron 
peak in the observed band, we have fitted routinely all datasets 
both with single and broken power-law models, with fixed 
Galactic and free absorption. The absorbing column density 
was parameterized in terms of N-n, the HI column density, 
with heavier elements fixed at solar abundances and cross sec- 
tions taken from Morrison & McCammon (1983). The Galactic 
value was derived from the nh program at HEASARC (based 
on Dickey & Lockman 1990), except for IBS 1218+304, for 
which a direct measure is available (Elvis et al. 1989). The re- 
sults from the single power-law fits are reported in Table 2, 
while Table 3 reports the broken power- law fits parameters. 
Adopting a conservative approach, we assessed the presence 
of a spectral break when the broken power-law model was sig- 
nificantly better, according to the F-test, even with respect to 



the single power-law model with free absorption. These are 
the cases presented in Table 3. Instead, when the fits were not 
significantly improved, we considered the single power-law 
model with free or galactic absorption as an adequate repre- 
sentation of the data. 

For most of the single power-law models, as shown in 
Table 2, the fits are improved (with significance > 99.9%) 
adopting an absorption value greater than the galactic one. In 
most cases, however, the broken power-law models with galac- 
tic absorbing columns provided even better fits (according to 
the F-test, see Table 3). For lES 0033+595 and PKS 0548- 
322, instead, there is evidence for a higher coluimi density also 
using the broken power-law models. To check if this extra ab- 
sorption is real, we have used for these datasets also the "curved 
spectrum" model, as in Fossati et al. 2000. 

In the following sections we present the details of the anal- 
ysis for each object, along with some notes on the individual 
characteristics. 
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ai 


-E'brcak 


a2 


Fi keV 


F2- 


10 kcV 




Lees 


Xr/d.o.f. 


F-test 


iu cm 




keV 




/iJy 


ergs cm 


s 


iviecs 






0033+595 


fits LECS+MECS+PDS 


















42.4 fixed 


0.34 ± 0.13 


2.4 ± 0.2 


1.03 ± 0.03 

i-U»_,_o.i2 


8.7 ±0.7 
13.2 ±2.5 


5.9 
6.0 


xlO" 
xlO" 


11 
11 


0.70 
0.70 


1.12/150 
1.05/149 


>68% 
>99.6% 


0120+340 


observation of 2/2/99 


















5.15 fixed 

^•"+4.9 






1.32 ±0.08 
1.32 ± 0.09 


*-0+1.9 


1.3 
1.3 


xlO" 
xlO" 


11 
11 


0.70 
0.70 


0.92/93 
0.93/92 


>97% 
>92% 


0548-322 


observation of 20/2/99 


















2.51 fixed 


3-°-^ 


1 o-°-2 

4 5-2-3 


1.04 ± 0.05 

1 4-0.3 


6 5-°-* 
5.7 ±0.3 


2.3 
2.3 


xlO" 
xlO" 


11 
11 


0.73 
0.77 


1.04/83 
0.95/82 


>97% 
>99.8% 


1114+203 


1.36 fixed 
1.5 ±0.6 


1 07-0.07 

^■■-"+0.06 

1 4n-°-22 

-^■^•^+0.18 


J^-"+0.4 

1 6"°-* 

■■■•"+0.7 


1.89 ±0.08 
1.89 ± 0.09 


4.4 ±0.4 
4.4 ± 0.4 


6.2 
6.2 


XlO" 
xlO" 


12 
12 


0.77 
0.77 


0.93/62 
0.95/61 


> 99.9% 

> 99.9% 


1218+304 


1.78 fixed 

± 0.7 


1 -ir-O.OS 
-'-•-'-^+0.07 

(J.l!) 


' +0.4 
i-'> 1 o.r, 


1.56 ± 0.05 

i. ."•)(; ± ().()() 


7 O""-'' 
' ■^'+0.6 

().!) ± ().() 


1.5 

i..") 


xlO" 

XlO 


11 

11 


0.72 

0.77 


0.84/92 

0.85/91 


> 99.9% 

> 90. 


2356-309 


fits LECS+MECS+PDS 


















1.33 fixed 
1.4 ±0.5 


0- 78+0.06 
O-80+o.ll 


-■-•8+0.4 

-■■o+o.s 


1.10 ±0.04 
1.10 ±0.05 


6.2 ±0.2 
6.2 ±0.2 


2.5 
2.5 


XlO" 
xlO" 


11 
11 


0.73 
0.73 


0.95/35 
0.97/34 


> 99.9% 

> 99.9% 



Errors at 68% conf. level for 3 and 4 parameters of interest (Ax^ = 3.52 and 4.72 respectively). 
Table 3. Broken power-law fits with galactic and free absorption. For these datasets, the broken power-law model is statistically 
preferred with respect to a single power-law model (the F-tests in the last column are against the single power-law model with 
free absorption). 



sorption). The best fit, however, still requires a column density 
in excess of the galactic value (see Table 3), with high signif- 
icance (F-test >99.8%). The spectral indices before and after 
the break energy are lower and higher than 1 , respectively, thus 
locating the synchrotron peak at the break energy itself, around 
3 keV. The PDS data agree well with the extrapolation of the 
power-law after the break. The obtained A^h values higher than 
the galactic one cannot be accounted for by a curved model: 
this gives unplausible values for the spectral index at low ener- 
gies (ai = 0.1 at 0.5 keV), for galactic absorption. More likely 
values (ai = 0.6 — 0.8) are obtained with a column density 
close to the one given by the broken power-law fit (5.7 — 6.3 
xlO^'^ cm"^). This object has been observed also by the Hubble 
Space Telescope as part of the "snapshot survey" of BL Lac 
objects, and was resolved into two point-Uke sources, with a 
separation of 1".58. It is a possible case of gravitational lens- 
ing (for a discussion of this hypothesis, see Scarpa et al. 1999, 
and Falomo & Kotilainen 1999). 

3.2.2. 1ES 0120+340 

This source was observed twice in a month, and was about 
^ 20% weaker during the second observation. The best fits are 
shown in Fig. 4. In both epochs, a good fit is obtained with the 
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Fig. 3. Broken power-law fit to the LECS+MECS+PDS data of 
lES 0033+595, with free galactic absorption. 



3.2.1. 1ES 0033+595 

This source was detected in the LECS instrument only above 
0.4 keV, due to the high galactic absorption, and in the PDS 
instrument up to ~ 60 keV. The spectrum is characterized by 
a curved shape, which is better fitted by a broken power-law 
than by a single power-law model (with both fixed and free ab- 
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Fig. 4. Best fits data and folded model, plus residuals, for the two observations of lES 0120+340. On the left, the model is a 
single power-law + galactic absorption (see text); on the right a broken power-law + galactic absorption. PDS data are probably 
contaminated (see text), and are not used for the fits. 
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Fig. 5. Confidence contours for the photon index and column 
density from the single power-law fits to lES 0120+340 data. 
The upper one is for the 2/2/99 observation, the lower for the 
3/1/99 one. The Une indicates the galactic column density. 



single power-law model, but with a free absorption well in ex- 
cess of the galactic value, at more than 99.9% confidence (see 
Fig. 5). As shown in Fig. 5 (and also indicated by the different 
ratio of MECS/LECS countrates), the spectral index has varied 
significantly, so it was not possible to sum the observations to 
enhance the S/N ratio. 

In the second observation, a better fit is provided by a bro- 
ken power-law model with galactic absorption, according to 
the F-test (>97%). This is also consistent with previous results 
by Brinkmann et al. 1994 (RASS data), which reported absorp- 
tion values near the galactic one [a = 0.93 ± 0.35, iVn = 
5.1±1.4x 10^0 cm-2, xl= 0.97(13d.o./.), = 4.50x10-1^ 
erg s~^ cm^^ in the ROSAT band]. In the first observation, in- 
stead, a higher column density is required even with a broken 



power-law model [x^=0. 98(120 d.o. f.), to be compared with 
X^=l. 10(121 d.o.f.) found with galactic Nh]. Comparing the 
two LECS datasets, the main differences lie in the 0.1-0.3 keV 
range, where oddly the higher flux observation presents a lack 
of soft photons with respect to the lower flux data. A check on 
the local background, compared to the blank-sky one, showed 
no relevant differences. Although this behavior could in princi- 
ple be due to a variable absorption, either local or at the source, 
the low S/N in that energy range does not allow to draw reliable 
conclusions. Beyond 0.3 keV, the galactic column density pro- 
vides good fits both with single and broken power-law models. 
Given the previous ROSAT results, we regard it as the most 
likely hypothesis, however future observations with higher S/N 
are necessary to settle this issue. 

The source was detected also with the PDS (see Fig. 4), 
but with a flux incompatible with that expected by a simple 
extrapolation of the LECS+MECS data. Possible explanations 
are a flattening (a ~ 1) of the spectrum due to the emerging 
of an inverse Compton component (not unconceivable, given 
its SED, see Fig. 14), or the contamination by other sources in 
the field of view not visible in the imaging instruments. This 
last hypothesis is corroborated by the presence, at 36' from the 
BL Lac object, of flie galaxy NGC 513 (z=0.0195), identified 
as a Seyfert 2. Objects of this class, being heavily absorbed, 
can actually have a significant hard X-ray emission, even if the 
soft X-ray flux is low (the ROSAT WGACAT fuU band coun- 
trate for NGC 513 is 0.0097 ± 0.0014 cts/s, corresponding to 
1.45 xlO~^'' erg cm^^ s^^ for galactic absorption). To check 
this possibility, we have also added a power-law component 
to the model for the PDS data of lES 0120+340, to mimic the 
Seyfert 2 emission in the hard band, after accounting for the 
off-axis efficiency of the instrument. In this way, the PDS data 
of both epochs can be easily explained by a power-law compo- 
nent of a ^ 0.7 with a normalization at 1 keV of 2.5 ± 1.4 /iJy, 
which requires a column density of « 4 x 10^^ cm~^ to account 
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Fig. 6. Best fit data and folded model, plus residuals, for the first and third observations of PKS 0548-322. The left panel shows 
the best fit model using a broken power-law model, the right panel using a single power-law, both with free absorption. 
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Fig. 7. Confidence contours for the photon index and column 

density from the single power-law fits to PKS 0548-322 data. 
The upper one is for the 7/4/99 observation, the lower for the 
20/2/99 one. The fine indicates the galactic colunm density. 



for the ROSAT countrate. These values are quite common for 
Seyfert 2 objects, so this is a likely possibiUty. 

3.2.3. PKS 0548-322 

The planned observation for this object has been split in three 
parts. The second observation (26-27/2/99) was very short (~ 
2000 s), and only MECS data are available. For the first ob- 
servation (20/2/99), an extraction region of 6' of radius has 
been used for the LEGS, due to the presence of a contaminat- 
ing source in the field of view, at ~ 10'. We have identified 
this source with the star GSC_07 61_01558 in the Guide 
Star Gatalog, probably in a flaring state: it is much less visible 
in the MEGS image (indicating a very soft spectrum) and was 
completely absent in the LEGS and MEGS images of the third 
observation. 



In both observations with LEGS data, a single power-law 
fit with free absorption gives values higher than the galactic 
one, at more than 99.9% confidence (see Fig. 7). The spectral 
index has varied significantly between the three observations, 
not allowing the sum of the datasets to enhance the S/N. For 
the third observation the single power-law model fits the data 
very well; instead, the first observation is significantly better 
fitted by a broken power-law (F-test > 99.8%), with a break 
around 4 keV. In both cases, the data require a column density 
in excess of the galactic one, with high significance (F-tests > 
99.9% and > 99% respectively, see Table 2 and 3). A colunm 
density higher than galactic seems also suggested by the fit with 
the "curved" model: the fit with a galactic Nh provides low 
energies spectral indices sUghtly flatter than expected from a 
peak around 3-4 keV (ai = 0.4 at 0.5 keV). More likely values 
are obtained with Nu around 3.6^.5 x 10^° cm"^ (ai = 0.7 - 
0.8). 

Interestingly, this extra absorption could be associated with 
the BL Lac environment. In fact, PKS 0548-322 belongs to a 
rich cluster of galaxies, and has two companions, one of which 
is seen in close interaction with the host galaxy of the BL Lac 
object (from NTT images, see Falomo et al. 1995). Previous 
X-ray observations with ROSAT and ASGA (see Sambruna 
et al. 1998) show some evidence of spectral absorption fea- 
tures around 0.5 keV, which suggest the presence of circum- 
nuclear ionized gas, perhaps fueled by the tidal interaction. 
Unfortunately our data, due to the low S/N ratio in the low 
energy band, cannot confirm nor exclude these features. 

PKS 0548-322 was detected up to ~ 30 keV by the PDS: 
the data are compatible with the extrapolation of the power-law 
from lower energies. 

3.2.4. GB 1114+203 

This source presents a very steep spectrum, with a further 
steepening beyond 1.6 keV best fitted by a broken power-law 
model. The column density is in agreement with the galactic 
value, within the errors. No detection was achieved in the PDS 
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Fig. 8. Best fit data and model for GB 1114+203: broken 
power-law with fixed galactic absorption. The source was not 
detected in the PDS. 
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Fig. 9. Best fit data and model for lES 1218+304: broken 
power-law with absorption fixed at the galactic values. 



instrument. Although the flux level at 1 keV is comparable with 
that of other sources in the sample, the steepness of the spec- 
trum {a 1.9 beyond 1.6 keV) makes the source too faint to 
be detected in the PDS range. The two steep indices locate the 
synchrotron peak at lower energies than the observed X-ray 
band, thus identifying this object as a typical "normal" HBL. 

3.2.5. 1ES 1218+304 

As the previous object, also this source presented a steep spec- 
trum best fitted by a broken power-law with high significance 
(> 99.9%). The found column density is in agreement with the 
galactic value. The source is detected in the PDS up to ^ 30 
keV, and the data, although slightly higher, are compatible with 
the extrapolation of the LECS+MECS best fit model within the 
(large) uncertainties (Fig. 9). Again, the two steep X-ray in- 
dices identify this object as a normal HBL. Previous EXOSAT 
and ROSAT observations reported similar spectral properties, 
with steep slopes near the BeppoSAX values (a « 1.5 ± 0.1 
and a = 1.22 ± 0.03 in the ME and PSPC instruments, respec- 
tively; see Sambruna et al. 1994, and Fossati et al. 1998). 

3.2.6. 1ES 1426+428 

The spectrum of this object is very well fitted by a single flat 
power-law model, with galactic absorption: there is no clear 
sign of bending in the shape of the residuals (see Fig. 10). 
Broken power-law models do not statistically improve the fit, 
and moreover give values of aiand a2both flatter than unity 
[they are also only marginally different (a2-ai~0.1)]. The 
data in the PDS band are particularly important for this source, 
in order to assess or not the presence of a steeper component 
(a > 1) that would reveal the position of the synchrotron peak. 

Unfortunately, the PDS data are contaminated by another 
known hard source, namely the quasar GB 1428+422 (Fabian 
et al. 1998). This object, characterized by a very flat spectral in- 
dex (a < 0.5), lies at 41' from lES 1426+428, thus inside the 
PDS f.o.v. but outside the MECS f.o.v.. For a lucky coincidence 



however, the quasar was observed by BeppoSAX on 4/2/99, 
four days before our observation. Thanks to the collaboration 
between the two proposing groups, we have been able to es- 
timate the relative inter-contamination. Furthermore, we have 
also checked in NED and WGACAT databases for other poten- 
tially contaminating sources within 80' from lES 1426+428, 
selecting hard sources (HR > 0.8) and/or with a high countrate, 
excluding those already identified as stars, and estimating their 
flux in the PDS range under the assumption of a power-law 
spectrum. 

Besides the quasar, the highest contributions came from 
two sources: WGA J1426.1+4247 (a ~ 0.57, ROSAT coun- 
fi-ate =0.011 cts/s) and CRSS 1429.7+4240 (a ~ 0.87, ROSAT 
countrate =0.015 cts/s). Their estimated high energy X-ray flux 
turned out to be roughly two orders of magnitude lower than 
lES 1426+428 in the PDS band, so we did not consider tiiem, 
leaving GB 1428+422 as the only contaminating source. 

The PDS data of both blazars were contaminated one by 
the other. The strategy to disentangle the two contributions 
has been the following. First we fitted only the LECS+MECS 
data of both sources. Then we extrapolated these best fit mod- 
els into the PDS range, for estimating the relative contami- 
nations. We fitted flie combined LECS+MECS+PDS data of 
one source, adding in the PDS range the LECS+MECS model 
of the other source, after accounting for the instrument off- 
axis response, and cross-checked the results for consistency. 
In Fig. 1 1 the single components of the models used in the lES 
1426+428 fit are shown, together with the total expected flux 
in the PDS instrument. Of course, this technique is based on 
the assumption that no large changes of flux or spectral index 
occurred in the four days separating the two observations, but 
in any case it gives a reference point to estimate possible vari- 
ations. Fig. 10 (left panel) shows the fit including for the PDS 
model flie BeppoSAX spectrum of GB 1428+422 (a = 0.42, 
Fi kcV = 0.30/iJy). The PDS points remain slightly above 
die flat (a^ = 0.92) power-law fliat fits flie lES 1426+428 
LECS-I-MECS data, however flie fit is good (xl = 1.06), witii 
no sign of declining slope up to 100 keV. A broken power-law 
model gives even flatter indices above 10 keV (a2 « 0.3 — 0.4), 
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Fig. 10. Single Power-law fits to the lES 1426+428 data. The PDS data have been deconvolved with the GB 1428+422 normal- 
ization either fixed at BeppoSAX value (left panel) or free to vary (right panel). In both cases the slope of the GB 1428+422 
spectrum is fixed at the BeppoSAX value {a = 0.42). 
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Fig. 11. The two models used to fit lES 1426+428 PDS data 
(shown in Fig. 10). Together with the model for the pointed 
source (the power-law with a = 0.92), the off-axis GB 
1428+422 contribution is shown, with the two different flux 
levels used. The lower one (thick fine) refers to the BeppoSAX 
observation value, the upper (dashed line) to the level resulting 
from the best fit obtained with the quasar flux as free param- 
eter. The two curved Unes correspond to the total flux in the 
PDS range, i.e. the sum of the a = 0.92 power-law with each 
of the two GB 1428+422 flux levels (thick and dashed lines 
respectively). 
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Fig. 12. Broken power-law fits to the LECS+MECS+PDS data 
of H 2356-309. Absorption fixed at the galactic value. 



flie PDS data letting the GB 1428+422 normalization free to 
vary (Fig. 10, right panel). Again, the lES 1426+428 spec- 
trum remains flatter than unity up to 100 keV, even account- 
ing for a GB 1428+422 increase of a factor ~4 in 4 days, from 
Fi kcv 0.30 /iJy to ^1 kcv = 1-44 ^Jy. Given the SED 
of lES 1426+428 (see Fig. 14), fliis result quite reliably con- 
strains the synchrotron peak to lie near or above 100 keV. 



however, these values are not very plausible, since, even not 
considering the GB 1428+426 contribution, the emerging of 
such a flat component would have imphed a 3 a detection in the 
PDS also in the 100-200 keV band, not observed (~ 8.2 x lO'^ 
cts/s against the observed 2.2 ± 2.7 x 10^^ cts/s). 

Assuming also the possibility of a large flux variation for 
GB 1428+422 (not unUkely: in 1998, during pointed ROSAT 
HRI observations, it has doubled its flux over a timescale of 
two weeks or less, see Fabian et al. 1999), we also fitted 



3.2.7. H 2356-309 

The spectrum of this source is not compatible with a single 
power-law model: there is clear evidence for a convex spectral 
shape, that is best fitted by a broken power-law with high sig- 
nificance (> 99.9%, see Table 3). This locates the synchrotron 
peak around 1.8 keV. The found column density is around the 
galactic values. This source has been detected also in the PDS 
instrument up to ^ 50 keV (see Fig. 12) and the data lie on the 
extrapolation of the MECS power-law. 
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4. Comparision with previous X-ray 
observations 

Among the seven sources of our sample, PKS 0548-322 and 
lES 1426+428 have been observed several times also by previ- 
uos X-ray missions. Here we compare our BeppoSAX results 
with older X-ray spectra, in order to reveal and study possible 
variations of the synchrotron peak frequency. 

4.1. PKS 0548-322 

Variations of the location of the synchrotron peak have been 
observed even during our BeppoSAX observations, showing a 
convex spectrum peaked around ~ 4 keV in the first dataset 
and a steep spectrum in the third observation (locating the syn- 
chrotron peak below the observed band). The ASCA observa- 
tions (Sambruna et al. 1998) confirm the presence of a con- 
vex spectrum, but with the break at slightly lower energies, 
around 1.7 keV (asoft = 0.9 ± 0.1 and ahard = 1-04 ± 0.02). 
The ROSAT slope agrees with the lower energy index of the 
ASCA broken power-law (a = 0.84 ± 0.05, Sambruna et 
al. 1998). The EXOSAT and GINGA satelhtes observed the 
source three and five times respectively, showing again an in- 
teresting spectral variability. The EXOSAT observations of 2 
November 1983 (GariUi & Maccagni 1990) and 7 March 1986 
(Barr et al. 1988) report a "convex" spectrum similar to our 
first dataset, with break energies indicating a synchrotron peak 
at 2.8 ± 0.8 and ~ 5 keV, respectively. The 30 Nov. 1983 ob- 
servation, instead, is well fitted by a steep single power-law 
(a ^ 1.2t"-2 in the 0.02-7.2 keV band), as our third dataset 
(GariUi & Maccagni 1990). Four out of the five GINGA obser- 
vations (Tashiro et al. 1995) carried out between January 1990 
and March 1991, showed an X-ray slope ^ 1 between 2-30 
keV (hence the synchrotron peak lies in this range). One obser- 
vation instead (19 Feb. 1991), quite interestingly, was charac- 
terized by a very flat spectral index (a = 0.84 ± 0.02) over the 
whole energy range, which would locate the synchrotron peak 
at energies >30 keV. These data suggest a possible behavior 
similar to Mkn 501, with the synchrotron peak which some- 
times reaches the hard X-ray band, even if most of the time is 
rather stable in the range 1-5 keV. 

PKS 0548-322 has been observed also in the TeV band, but 
so far no detection was reported: CANGAROO observations 
between October and November 1997 showed no evidence for 
an on-source excess (Roberts et al. 1999). The simultaneous 
optical and XTE All Sky Monitor X-ray measurements did not 
show any special activity. 

4.2. 1ES 1426+428 

Besides ROSAT, this source has been previously observed also 
by HEAO-1, EXOSAT, ASCA and the BBXRT experiment 
(details in Sambruna et al. 1994, 1997). In Fig. 13 we report 
the historical flux levels in the 2-10 keV band, together with the 
measured spectral indices in that band. In all the previous ob- 
servations, there is no evidence of a flat spectrum above 2 keV, 
even in the presence of large flux variations. The BeppoSAX 
observation is the first where this source showed a flat spec- 
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Fig. 13. Historical X-ray fluxes (2-10 keV) measured by the 
satelhtes HEAO-1, EXOSAT, ROSAT, ASCA, and flie BBXRT 
experiment (see Sambruna et al. 1997). The fluxes have been 
obtained from single power-law fits with galactic absorption 
to the 2-10 keV continuum (except ROSAT, where the values 
are from the extrapolation in the 2-10 keV band of the best 
fit model). In the upper part are also displayed the correspond- 
ing spectral indices. The errors are estimated from those on 
normalization and spectral index, at 90% confidence level. For 
ROSAT, ASCA and BeppoSAX , they are comparable with the 
marker size. Note how in the BeppoSAX observation the flux 
level was similar to the ASCA and ROSAT ones, but with a 
very different spectral index, indicating a shift in the peak fre- 
quency of about two orders of magnitude. 



trum up to 100 keV. In BBXRT and ASCA data, the full band 
spectrum is best fitted by a broken power-law: in this case, the 
resulting break energy locates the synchrotron peak, for both 
observations, between 1 and 2 keV. This source therefore has 
undergone a shift in the peak frequency of at least two orders of 
magnitude. Note that, as shown by Fig. 13, the BeppoSAX flux 
in the 2-10 keV band of 1 ES 1 426+428 is not particulariy high: 
it is the second lowest state observed, about ~20% higher than 
the lowest one, observed by ASCA (2.04 vs. 1.64 x 10~^^ ergs 
cm^^ s^^, in the 2-10 keV band). This indicates that the object 
can be very powerful above 10 keV while remaining inconspic- 
uous in the soft X-ray band. This is also supported by the XTE 
monitoring, which did not show any special activity during the 
BeppoSAX observation. The other two "over 100 keV" sources 
(i.e. Mkn 501 and lES 2344+5 14) showed a markedly different 
behavior, becoming very active and strong also below 10 keV 
when their synchrotron spectrum peaks at high energies. lES 
1426+428 has shown that the soft X-ray activity is not always 
necessary to be a strong 100 keV synchrotron emitter. At the 
same time, in this source, a relatively faint 2-10 keV state does 
not guarantee a hard spectrum (see e.g. the ASCA and ROSAT 
data in Fig. 13): all this is revealing of a complex relation be- 
tween the injected power and the location of the synchrotron 
peak. 



L. Costamante et al.: Extreme Synchrotron BL Lac Objects 



11 



I 

m 
I 



O 
Cifl 

u 



o 



-10 

-11 

-12 
-13 
-14 

-11 
-12 
-13 
-14 
-15 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



o 



0120+340 
z = 0.272 



ID 



n I I I I I I I I I I I I I I 



o -- 



1426+428 
z = 0.129 



_I I I I I I I I L 



-| — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

o 



o 



o 



0548-322 
z = 0.069 



I I I I I I I I I I I I I I I 



:: o 



2356-309 
z = 0.165 

J I I I i_ 



10 15 20 25 10 15 20 25 

Log V [Hz] Log V [Hz] 



Fig. 14. The SEDs of the 4 BL Lacs observed in A02, made with BeppoSAX (small filled points) and literature data (open circles, 
bow-ties for slope informations and arrows for upper limits). All four sources peak in the X-ray band. 



5. Discussion 

The main purpose of our observations was to characterize the 
global spectral properties of the selected candidates, verifying 
their "extremeness" by locating the position of the synchrotron 
peak by means of the X-ray spectrum itself and by the shape of 
the overall SED. The SEDs of the 7 sources, constructed with 
our BeppoSAX and literature data, are shown in Fig. 14 and 
15. The BeppoSAX observations have revealed and confirmed 
the "extreme" nature of five objects: for four of them (lES 
0033+595, lES 0120+340, PKS 0548-322 and H 2356-309) 
the X-ray data themselves show the peak of the synchrotron 
emission, presenting a convex spectrum best fitted with a bro- 
ken power-law whose break locates the peak in the 1-5 keV 
range. For lES 1426+428, instead, the flat X-ray spectrum up 
to 100 keV constrains the synchrotron peak to lie near or even 
above that value, and makes this source the third object ever 
found with such high peak frequencies, after Mkn 501 and lES 
2344+514. For the last two sources (GB 1114+203 and lES 
1218+304), again best fitted by a broken power-law, the in- 
dices are both steep, thus locating the peak below the observed 
X-ray band (see Fig. 15). Their steep spectrum qualifies them 
as typical HBLs, for which the X-ray band is dominated by the 
tail of the synchrotron emission after the peak. These results 
confirm that the criteria adopted to select our candidates were 
highly efficient (5 out of 7), and we are beginning to populate 
the high energy end of the synchrotron peak sequence. 



An important aspect that our observations have pointed out, 
compared with older data, is the high variability of the peak 
frequency exhibited by these objects: the observed range of 
j^pcak variations has been from less than 0.1 to 1.4 keV for lES 
0120+340, from < 0.1 to > 20 keV for PKS 0548-322, and 
from -1 to over 100 keV for lES 1426+428. This can affect 
the classification of an object as "extreme": in fact, if observed 
only once by BeppoSAX , lES 0120+340 and PKS 0548-322 
would have been classified as normal HBLs (see Table 2). Mkn 
501 and lES 2344+514 too exhibited HBL properties before 
their observations during flaring states. It is likely, therefore, 
that other "extreme" objects could be hidden in the HBL class, 
and would be uncovered if observed many times. Repeated ob- 
servations are therefore necessary to properly identify this type 
of objects, and also to sample their spectral behavior in connec- 
tion with flux changes. Flux variations, however, should not be 
the only criteria to trigger the observations for the most ex- 
treme states: lES 1426+428 proves that the shift in the peak 
frequency can be de-coupled from high flux states in the soft 
X-ray band. 

A modeling of the SEDs and the resulting physical parame- 
ters will be presented in a forthcoming paper (for a preliminary 
Synchrotron Self Compton model to the A02 sources, see Fig. 
1 of Costamante et. al. 1999). Although more objects are neces- 
sary to draw reliable and statistically significant conclusions, it 
is nevertheless interesting to start analyzing the characteristics 



12 



L. Costamante et al.: Extreme Synchrotron BL Lac Objects 





-10 


' 1 ' 


1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


' 1 ' 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 








vj 


(§> ^ - 






- 


1 

m 


-12 




O 










1 
















b 


- 14 






0033+595 

\J \J K-f 1 kJ l_/ 






1218+304 


r \ 






Z — 


- n 




Z — U . 1 OiC 


U 




H — 1 — 1 — 1 — 


1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — \ 


"9i , 






-10 


— 1 — 1 — 1- 






10 


15 


20 25 








o 


i-©i 

\ 




Log y 


[HzJ 




-12 














:\ 


















-14 






1114+203 J 








o 


, 1 , 


, , , 1 


z= ? ; 

r . . 1 . . . . 1 . 1 












10 


15 


20 25 














Log 


y [Hz] 









Fig. 15. The SEDs of the 3 BL Lacs observed in A03, made with BeppoSAX. and hterature data. Only IBS 0033+595 peaks in 
the observed X-ray band; the other two objects peak at lower energies. A polynomial fit as in Wolter et al. 1998 ( \og{vFy)=a + 
h{\ogv)+c{\ogvf + d{\ogvf) gives log(i/peafc)=15.4 and 16.3 for GBl 1 14+203 and lES 1218+304, respectively. 



of these extreme objects also in the context of the entire class 
of BL Lac objects. 

Fig. 16 shows the broad band spectral index a^x as a func- 
tion of the frequency of the synchrotron peak, for the BL Lacs 
analyzed in this paper plus Mkn 501 and IBS 2344+514, to- 
gether with the HBLs analyzed in Wolter et al. (1998) and the 
RBLs analyzed in Comastri et al. 1995. It can be seen that 
the almost linear correlation between arx and log v-pcak fat- 
tens out for large values of i^pcak, reaching a "stable" value of 
Urx ^ 0.5. This is in agreement with the scenario of a syn- 
chrotron peak moving smoothly from lower to higher energies: 
when the peak moves beyond the X-ray band, both the radio 
and the X-ray fluxes are produced by the same branch of the 
synchrotron emission, and so both fluxes change similarly as 
the peak moves at still higher frequencies. Consequently a^x 
stabilizes at a (flat) value. Instead, when the radio and X-ray 
fluxes are produced by different branches of the synchrotron 
emission (before and after the peak), the radio and the X-ray 
fluxes change differently as the peak shifts, thus changing arx- 

The "moving peak scenario" also nicely accounts for the 
relation between ax and t'pcak (Padovani & Giommi 1996), 
shown in Fig. 17. To reliably compare our results with pre- 
vious ROSAT results, in Fig. 17 we plot BeppoSAX spectral 
indices in the soft band (i.e. the first spectral index for a bro- 
ken power-law model). We can see that LBLs are characterized 
by flat indices (ax < 1), since their X — ^ray emission is likely 
dominated by the flat inverse Compton process. Intermediate 



BL Lacs, with I'poak ^ 10^''-10^^ Hz, are characterized by 
Qfa; ~ 1. For them, broad band X-ray observations should re- 
veal a concave spectrum, with the soft X-rays dominated by 
the steep tail of the synchrotron emission, and with the flat in- 
verse Compton radiation emerging at higher energies (see the 
nice example of ON 231 in Tagliaferri et al. 2000). The steep- 
est indices correspond to HBLs, where the X-rays are entirely 
dominated by the steep tail of the synchrotron emission. But 
as t^pcak increases further, the synchrotron peak moves into the 
X-ray band, making the spectrum to flatten again, eventually 
reaching a "stable" value (ax ^ 0.5) corresponding to the typ- 
ical slope of the synchrotron spectrum well before its peak. 
These are the objects that we have called extreme BL Lacs in 
this paper (i.e. with i^pcak > IkeV). 

The left panel of Fig. 18 shows z^pcak as a function of the 
radio luminosity (vL^ at 5 GHz) for our targets and for a few 
samples of blazar objects (IJy BL Lacs, Slew Survey BL Lacs 
and 2Jy FSRQs), including the objects studied by Wolter et 
al. (1999) as weU as Mkn 501 and IBS 2344+514 in the flar- 
ing state. The solid line is the phenomenological relation pro- 
posed by Fossati et al. (1998) for the analytic representation 
of the general SBD. In the context of the samples under con- 
siderations, our extreme BL Lacs, having luminosities similar 
to HBLs but larger J^pcak, indicate a wider range of peak fre- 
quencies covered by HBLs, compared to lower peaked objects. 
Indeed, the synchrotron peak frequency of the extreme BL Lacs 
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Fig. 16. The radio to X-ray broad band spectral index arx as 
a function of i^pcak- Open symbols correspond to data from 
Comastri et al. 1995 (crosses and squares), and Wolter et al. 
1998 (circles). Filled circles are the sources analyzed in this 
paper The stars represent the quiescent and flaring states (con- 
nected by a Hne) for lES 2344+514 and Mkn 501. The arx 
parameters have been calculated as in Wolter et al. (1998), as 
well as the i^peak, when not directly measured in the X-ray 
spectrum. 

and some HBLs is highly variable (as in the case of Mkn 501, 
lES 2344+514 and lES 1426+428, see Fig. 18). 

At present, we do not know if this extreme spectral vari- 
ability is a characteristic common to all blazars or if it is a pe- 
culiarity of only the HBLs (or some of them). It is however 
interesting to note that for some well studied LBL and FSRQ 
the value of z^pcak seems to remain more constant. As an il- 
lustration we show in Fig. 18 (right panel) the observed range 
of variation of both j/pcak and Lpoak = i^poak-^i^pe^k for three 
"low peak" sources of different luminosity (BL Lac itself, OJ 
287 and 3C 279). The high //peak objects seem to be character- 
ized by a much larger i^poak variability, compared to the others, 
given more or less the same range of luminosity variations. 

The high synchrotron peak frequencies which characterize, 
by definition, the extreme BL Lac class, flagging the presence 
of high energy electrons, make these sources good candidates 
for Te V emission through the inverse Compton process, as also 
indicated by the TeV detections of Mkn 421, Mkn 501 and 
lES 2344+514 (Catanese et al. 1997, 1998; see also Stecker, de 
Jager & Salamon 1996). A significant role in the propagation of 
such high energy photons, however, is played by the absorption 
resulting from pair production interactions with the intergalac- 
tic infrared background. In fact, it has been shown that photons 
with energies greater than ^ 1 TeV should be removed by the 
spectra of sources with redshifts > 0.1, according to the esti- 
mates of the intergalactic IR background (Stecker & De Jager 
1997, 1998). In this respect, the extreme BL Lacs could be very 
useful as probes for the IR background, thanks to their different 
redshifts which make possible a good sampling of the optical 
dephts along z. IBS 1426+428, in particular, should be very in- 
teresting, given its high synchrotron frequencies (in the "Mkn 




Fig. 17. The X-ray spectral index ax as a function of t'poak 
(Padovani & Giommi 1996). Open squares and circles are from 
Comastri et al. (1995) and Wolter et al. (1998). The values from 
Comastri et al. (1995) are obtained from ROSAT observations 
(for BeppoSAX results on LBLs, see Padovani et al. 2001). 
Filled circles mark the BL Lacs analyzed in this paper For con- 
sistency, the ax here reported is the soft one, in case of broken 
power-law best fits. The stars (Mkn 501) and asterisks (IBS 
2344+514) refer to the quiescent and flaring states values for 
these two sources. 

501 class") and its redshift at the border of the expected ob- 
servable region (z — 0.129). 
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